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Введение. При организации бурильных работ одна из основных проблем — точность и плавность опускания 
связок труб в шахту бурильной установки. Это зависит от многих факторов, в том числе от работы 
гидротормоза подъемного устройства. Цели данной работы: создание и исследование математической модели 
гидродинамического балласта в буровой установке. С помощью метода обратного анализа изучено влияние 
некоторых эксплуатационных показателей на тормозной момент гидравлического тормоза. 

Материалы и методы. Эксперименты проводились с использованием лабораторной установки, которая 
представляет собой модель гидродинамического тормоза. Его клапан закрывали при различных условиях, 
чтобы получить несколько значений давления с расчетом тормозного момента при подвешенном грузе 
определенной массы. Моделировались реальные (полевые) условия работы гидродинамического тормоза и 
сравнивались полученные результаты. При создании математической модели использован метод обратного 
анализа. Он основан на результатах экспериментальных измерений и позволяет выразить совокупность влияний 
отдельных переменных на тормозной момент. 

Результаты исследования. Создана и протестирована математическая модель гидравлического тормоза. 
Определена зависимость тормозного момента от давления, плотности и вязкости балластной жидкости. 
Влияние каждой переменной определено экспериментально, т. к. рассматриваемую зависимость нельзя 
представить в виде прямой связи. Методом обратного анализа получен набор значений констант, которые дают 
оптимальное решение. Принимая во внимание стандартный массив ошибок и минимальную стандартную 
ошибку, учтены статистические ошибки, допущенные во время экспериментальных измерений. 
Визуализирован физически приемлемый диапазон значений предложенной математической модели. С 
помощью базовой (нелинейной) математической модели рассчитан вспомогательный тормозной момент 
гидродинамического тормоза в зависимости от давления, плотности и вязкости. Установлена адекватность 
предложенной модели. В качестве критерия корректности использовались расчетные показатели тормозного 
момента. Ошибочное расхождение не превысило 6 %. Для дополнительного тестирования модели проведен 
расчетный эксперимент, имитирующий полевые условия. 

Обсуждение и заключения. Для математического моделирования гидродинамического балласта в буровой 
установке целесообразно использовать метод обратного анализа. Модель, предложенная в данной работе, 
связывает тормозной момент гидродинамического тормоза с рабочими параметрами жидкости внутри балласта: 
давление, вязкость и плотность. Доказана объективность модели. Предложена поправка к ней для симуляции 
работы тормоза в полевых условиях. Опираясь на полученные результаты, в будущих исследованиях 
целесообразно испытать созданную модель в полевых условиях с реальной полезной нагрузкой. 
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Введение. При обустройстве буровых установок операции подъема и опускания бурильных труб 
обеспечиваются подъемными устройствами (катушками). Они снабжены тормозами двух типов, основной из 
которых — механический. В процессе бурения вес труб превышает полезную нагрузку. Для его компенсации 
вводят вспомогательный регулировочный балласт, предназначенный для снижения скорости при установке 
группы бурильных труб [1-4]. 

Балласты используют для контроля нагрузки и скорости движения крюка, а также для поглощения 
кинетической энергии группы скважинных бурильных труб. Помимо этого балласты: 

— снижают усилия на буровой установке, особенно при больших нагрузках; 

— уменьшают износ основных механических балластных элементов; 

— помогают медленно и плавно остановить груз, прикрепленный к крюку. 

Один из типов вспомогательных тормозов — гидродинамический. В этом случае находящаяся в 
балласте вода преобразует часть поглощенной кинетической энергии подъемной оси в тепловую, и за счет 
этого опускаются и поднимаются трубы [5-7]. 

При манипуляциях с трубами ось тормоза соединяется с осью подъемных устройств. Во время 
операции подвижная часть установки вытесняет воду в направлении наклонных лопастей внутри статора и 
вращается со скоростью, равной скорости вращения оси подъемных устройств. Вода попадает на ребра статора, 
а затем на его детали. Процесс повторяется, и увеличиваются силы, препятствующие движению ротора. Это 
создает тормозной момент, который отражает вращательное движение осей подъемных устройств и уменьшает 
скорость вращения. В результате скорость спуска группы бурильных труб уменьшается [8-10]. 


'Антибас И. Р. и др. Метод обратного анализа при математическом моделировании гидродинамического балласта 





Коэффициент безопасности торможения (если его значение не превышает допустимого) координирует 
вспомогательную гидродинамическую работу и эффекты от работы двух основных тормозных систем, а также 
обеспечивает более длительное время обслуживания элементов главного тормоза. Это происходит за счет 
уменьшения износа фрикционных дисков и плоскостей тормозных колес. Увеличение тормозного момента 
балласта позволяет обеспечить правильное тормозное движение цилиндрического подъемника. Это 
свидетельствует о важности исследования гидродинамического тормоза. 

Гидравлическая мощность тормоза М определяется из соотношения: 

№=р.8-О.Н, (0 
где р — плотность рабочей жидкости (пресной воды); > — ускорение свободного падения; О — количество 
израсходованной рабочей жидкости, равное объему рабочей жидкости, которая проходит через систему ребер 
за один рабочий цикл; Н — высота (уровень) рабочей жидкости в балласте. 

Тормозной момент М» принудительного торможения двигателя определяется соотношением 

М, =р-8-О-Н/®, (2) 
где « — угловая скорость вращения движущегося участка балласта. 


Показатели кинематического торможения оребренных гидравлических машин определяются по 
известным теоретическим соотношениям тормозного момента: 


М, =^„-р-(5*-4°)-в, (3) 
и? 
М, = и -р-(`-4) о. (4) 
Здесь р — наружный диаметр «кольца» рабочей жидкости, образованного при вращении колеса ротора и 


принятого равным диаметру колеса подвижной секции; 4 — внутренний диаметр кольца рабочей жидкости, 
который зависит от уровня балласта; и — число оборотов движущегося участка балласта (число вращающихся 
осей), об/мин: 
2пи т 
® = > = 
60 100 

Коэффициент гидравлического тормозного момента Ам — безразмерная величина, учитывающая форму 
рабочих полостей в балласте, параметры и количество тормозных ребер. В практических расчетах среднее 
значение коэффициента гидравлического тормозного момента принимается 0,3. 

Цель исследования — изучить влияние некоторых эксплуатационных показателей на тормозной 
момент гидравлического тормоза, который представляет собой давление внутри тормозной камеры, а также на 
плотность и вязкость рабочей жидкости внутри балласта. Этих показателей нет в соотношении (4), поэтому 
примем, что тормоз работает при атмосферном давлении в 1 атмосферу, жидкость представляет собой пресную 








воду с плотностью 1 г/см” и вязкостью 1 стокс. 

Материалы и методы. Лабораторные эксперименты проводились в университете Алеппо (Сирийская 
Арабская Республика) на установке, представляющей собой модель ингибирования гидродинамического 
тормоза (рис. 1). 





Рис. 1. Лабораторная установка — модель ингибирования гидродинамического тормоза: | — шкив, 2 — подлокотник 
шкива, 3 — манометр, 4 — выпускной клапан, 5 — термометр, 6 — отверстие в верхней части бака, 7 — регулирующие 
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клапаны уровня жидкости, 8 — выходной канал для жидкости, 9 — выход жидкости из резервуара, 10 — клапан слива 
жидкости, 11 — впускной клапан, 12 — входная линия, 13 — гидравлический тормоз, 14 — катушка, 15 — измеритель 
количества циклов, 16 — трос, 17 — полезная нагрузка, 18 — концентратор связи, 19 — линия управления, 20 — бак для 
жидкости 
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Этапы экспериментов 

1. Резервуар и балласт заполняют жидкостью до требуемого уровня (0,106 м до уровня первого 
клапана). 

2. С помощью ручного рычага, связанного с катушкой, на 0,317 м поднимают подвешенный на крюке 
груз. Он весит 8 кг и соединен с тросом. 

3. Позволяют связке труб упасть под собственным весом. 

4. Снимают показания с индикатора оборотов оси катушки. 

5. Затем фиксируют влияние показателей, полученных в ходе лабораторных экспериментов по 
изучению гидравлического торможения: 

— давление, создаваемое в тормозной камере (Р); 
— плотность рабочей жидкости (р); 
— вязкость рабочей жидкости (1). 

Сначала рассчитывают тормозной момент вспомогательного гидродинамического тормоза по 
соотношению (4). С этой целью определяются значения и, Л, 4. 

Наружный диаметр «кольца» рабочей жидкости О образуется при вращении колеса подвижной секции 
(ротора) и принимается равным его диаметру. Диаметр Р движителя лабораторной установки в практических 
экспериментах равен 0,33 м. 

Внутренний диаметр кольца рабочей жидкости 4 зависит от уровня балласта. Для определения этого 
диаметра выравниваются объемы: 

— рабочей жидкости внутри балласта; 
— жидкости в «кольце», которое образуется при вращении подвижной секции внутри балласта и 
ограничивается высотой охлаждающей жидкости в баке (0,106 м). 

При Д = 0,33 м рабочий объем балластной жидкости — 2,85 литра. Расчеты показали, что внутренний 
диаметр 4 = 0,283 м. 

Число оборотов и связано с подвешенным грузом. Если тормоз не работает, то это число оборотов 
намоточного устройства. Если же оно отделено от тормоза, то и берется из лабораторных экспериментов: при 
весе груза 8 кги = 240 об/мин. 

Рассчитаем момент гидравлического торможения для веса 8 кг, прикрепленного к крюку. Во время 
эксперимента бак был заполнен пресной водой до уровня первого клапана (0,106 м) при следующих 
показателях: давление Р = 1 атм, плотность жидкости р = 1 г/см? и вязкость жидкости и = 1 Ст. Поэтому момент 
гидравлического торможения: 

2 
М, =0,3-1000- (0,33° 59,288). 29 =355 Нм. 
100 

Рассмотрим соотношение, обратное (2). В первую очередь речь идет об обратной зависимости между 
крутящим моментом и угловой скоростью (©). Уменьшение числа оборотов и на 20 % при увеличении 
плотности жидкости означает увеличение крутящего момента на 20%. Учтем допустимое значение тормозного 
момента 355 Н`м. В итоге получим искомую величину тормозного момента при следующих показателях: 
давление Р = 1 атм, плотность р = 1 г/см? и вязкость жидкости и = 1 Ст. 

Применительно к данной работе отметим следующее. При определении тормозного момента 
гидравлического тормоза учитывают число оборотов катушки, связанной с осью тормоза (то есть во время 
торможения): уменьшение числа оборотов намотки на определенный процент означает такое же увеличение 
крутящего момента. 

Влияние давления внутри балласта на тормозной момент гидравлического тормоза. Внутри 
гидравлического тормоза создали давление жидкости, закрыв клапан на линии выхода жидкости. Провели 
лабораторный эксперимент и расчет тормозного момента при подвешенном грузе массой 8 кг. Результаты 
представлены в таблице 1. 

Таблица 1 
Изменение тормозного момента гидравлического тормоза при изменении давления внутри балласта 














Давление внутри балласта, атм 1 1,25 1,4 1,55 1,7 1,85 
Левое граничное значение крутящего момента, Н-м 105 104 103 100 95 90 
Правое граничное значение крутящего момента, Н`м 355 359 363 372 388 405 























На рис. 2 показано изменение тормозного момента при изменении давления внутри балласта. 
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Рис. 2. Изменение тормозного момента при разном давлении внутри балласта 


Влияние плотности рабочей жидкости на тормозной момент. В экспериментах использованы 
химические вещества, которые увеличивают плотность воды, не подвергая коррозии элементы лабораторного 
оборудования. 

Рассмотрены различные концентрации растворенных в воде веществ для получения балластной 
жидкости разной плотности. Условия те же: исследуется тормозной момент при подвешенном грузе 


8 кг (таблица 2). 























Таблица 2 
Зависимость тормозного момента от плотности жидкости, числа витков катушки 
и добавок в балластную жидкость 
Плотность жидкости 1 1,065 1,09 1,11 1,13 1,15 
Число витков катушки 105 104 102 99 95 91 
Соли хлорида натрия 
Тормозной момент 355 359 366 375 388 402 
Пищевой сахар 
Тормозной момент 355 359 363 368 378 388 


























Отмечено очень слабое изменение вязкости при увеличении концентрации соли хлорида натрия. 
Можно считать, ЧТо вязкость жидкости приблизительно равна 1 Ст. 
На рис. 2 показана зависимость тормозного момента от плотности рабочей жидкости при 


использовании соли хлорида натрия и пищевого сахара. 
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6) 
Рис. 3. Зависимость тормозного момента от плотности рабочей жидкости с введением 
соли хлорида натрия (4) и пищевого сахара (6) 


Влияние вязкости рабочей жидкости на торможение гидросистемы. Использовались различные 
комбинации концентраций растворенного в воде глицерина и силиката для получения балластных жидкостей, 
отличающихся по вязкости. Условия те же: исследуется тормозной момент при подвешенном грузе 
8 кг (таблица 3). 























Таблица 3 
Зависимость тормозного момента от вязкости жидкости и добавок в балластную жидкость 
Вязкость жидкости 1 16 22 29 
Число витков катушки 105 102 98 94 
Глицерин 
Тормозной момент 355 365 378 392 
Силикат натрия 
Тормозной момент 355 372 388 408 

















На рис. 4 показаны изменения тормозного момента с изменением вязкости рабочей жидкости при 
использовании глицерина и силиката натрия. 
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6) 
Рис. 4. Зависимость тормозного момента от вязкости рабочей жидкости внутри 
балласта с введением глицерина (а) и силиката натрия (6) 


Полученную таким образом почти линейную зависимость нужно проверить с помощью 
математической модели. 

В ходе исследования индекса вязкости установлено, что параметры плотности и вязкости мало зависят 
друг от друга и плотность не равна 1 г/см? (таблица 4). 



































Таблица 4 
Соотношение плотности и вязкости жидкости 
Вязкость жидкости 16 22 29 23 30 38 
Плотность 1,008 1,014 1,021 1,017 1,025 1,034 

Коэффициент влияния плотности не превышал 6 %. Его рассчитывали как разницу между плотностями 

вязкой жидкости и пресной воды. Например, если вязкость была 30 Ст, то плотность — 1,025 г/см”. 
Следовательно, коэффициент воздействия плотности: 

1,025-—1 


В 055 


Результаты исследования 

Создание математической модели гидравлического тормоза. Зависимость тормозного момента от 
давления, плотности и вязкости нельзя представить как прямую связь. Влияние каждой переменной на 
тормозной момент определяется экспериментально (рис. 2—4). 
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Метод обратного анализа. В научной и справочной литературе нет зафиксированных показателей 
взаимосвязи тормозного момента с давлением, плотностью и вязкостью. Поэтому был задействован метод 
обратного анализа [4]. Он эффективен для создания математических моделей на основе экспериментальных 
измерений. С его помощью можно показать, как совокупность названных выше переменных влияет на 
тормозной момент. Это можно выразить соотношением: 


ВЕЕР). (5) 


где Е— функция, которая связывает изучаемое явление Рс набором значений переменных и набором значений 
для констант математической модели явления. 

Обратный анализ позволяет найти множество значений предполагаемых констант модели С через 
обратную зависимость: 


РеР"( В.С), (6) 


т? 
где Р„ — множество экспериментальных значений для изучаемого явления [4, 5]. 
Считается, что прямое решение для отношения (6) невозможно. Поэтому следует задействовать 
итерационную систему, чтобы получить набор значений констант, которые дают оптимальное решение. В 
обобщении указывается набор значений Р для соотношений: 


ВЕР.) (7) 
|р.-2,|<а, (8) 
где 2: = набор рассчитанных значений для изучаемого явления; & — требуемая точность в соответствии с 


расчетом изучаемого явления. 
Предложенный метод учитывает статистические ошибки, допущенные во время экспериментальных 

измерений, ориентируясь на стандартный массив ошибок (5,) и минимальную стандартную ошибку. Таким 

образом, согласно статистическому распределению Гаусса зависимость плотности статистических данных: 


ее соток] 5 [(Р. -р,) С(р.-Б, )] | (9) 
Здесь С, — стандартный массив ошибок. Символ { означает матрицу: 
аа о 
В белы (ана о |4, (10) 
иной 0, ,›52 (и) 


где и — количество контрольных точек. 
р. является функцией Р, поэтому проблема связана с определением такого значения Р, чтобы О. дало 
максимальное значение Р.. 


С —1 
5, =(Р.-2,) "С: (р,-р,). (11) 

Минимальное значение Ку можно Получить при более чем одном наборе значений Р. Поэтому 
правильный набор параметров формируется с помощью фактического физического моделирования значения 


этих параметров (рис. 5) [4, 5]. 





Рис. 5. Определение физически приемлемого диапазона значений (зеленая область): 
5, — физически приемлемое математическое решение, 55, 53 — математические решения 
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Поэтому мы используем набор первичных значений для параметров Ро, которые выводятся из средней 
нагрузки физического поля с учетом матрицы стандартных отклонений 5. 


При использовании нового распределения Гаусса получаем следующий статистический параметр 


интенсивности: 
1 в уд: 
РЕВ воть -5[(Р-Р) с; Р-в]] (12) 
Здесь С, —Щ_ матрица стандартного отклонения: 
аи анод 0 
С. =| 0... Вы о |. (13) 
аи ная 52 (>) 


Проблема, выраженная в отношении (6), решается путем обратного анализа. При этом набор 
значений Р находится в пределах, обозначенных в (7) и (8). Таким образом, можно определить общую область 


параметров Р! и Р› по распределению Гаусса: 
ТРЕ р.Р = соп$1.ехр(-5), (14) 


т а Е 
5=5(Р.-5.) -Сл(р.-0,)+(Р-В)-С»(Р-В)|. (15) 
Итак, для решения нужно найти максимальные или минимальные значения 5 функции Г (Р). 
Единственный способ найти минимальное значение 5 — использовать цифровые методы, например 
Ньютона — Гаусса. Этот подход основан на преобразовании аналитических отношений в цифровые 
итеративные с учетом ошибки, допущенной из-за игнорирования некоторых ограничений в аналитических 


отношениях (рис. 6). 











Метод обратного анализа 
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Рис. 6. Алгоритм метода обратного анализа 


На следующем этапе работы с помощью базовой (нелинейной) математической модели вычислили 
вспомогательный тормозной момент гидродинамического тормоза. Он зависит от давления, плотности и 
вязкости: 


еее 6) 
М, =а-р"+Ь- р" +а-и+С. (17) 
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Математическая модель (17) является общей нелинейной моделью, если значения п1, и, из отличны от 
единицы. Множество констант для поиска — это а, Б, 4, а также пи, и›, из. 
После применения итеративного цифрового метода с точностью, равной 0,001, найдены параметры, 


представленные в таблице 5. 











Таблица 5 
Значения параметров предлагаемой модели 
Параметры п п пз а Ь а С 
Величина 1 1 1 58,5 284,6 1,29 10,62 


























Согласно таблице 5, исследуемое физическое явление можно представить как линейную модель (пт, и>, 
пз) =1: 

М, =58,5Р+ 284, 6р-+1,29и+10,62. (18) 

Математическая модель основана на результатах лабораторных экспериментов и подразумевает 
гомогенизацию единиц измерения в соответствии со значением констант. 

Адекватность математической модели. Первый этап определения адекватности модели — это 
испытание на торможение при: 

— давлении Р = 1 атм, 
— плотности жидкости р = 1 г/см*, 
— вязкости жидкости и = 1 Ст. 

В данном случае торможение было равно начальному тормозному моменту Мь = 355 Н`м. Это логично, 
учитывая экспериментальные значения. 

Второй этап: три лабораторных эксперимента с подвешенным на крюке 8-килограммовым грузом. Они 
кратко описаны ниже. 

Первый. Взяли максимальные значения переменных и тормозного момента, а затем выполнили расчет с 
использованием математической модели по соотношению (18). Определили коэффициент ошибок, сравнив 
экспериментальный и математический результат. 

Второй. Взяли случайные значения переменных, относительно далекие от экспериментальных 
значений и от расчетных для тормозного момента (по математической модели). Вычислили коэффициент 
ошибок, сравнив экспериментальный и математический результат. 

Третий. Взяли случайные значения для переменных, относительно далекие от экспериментальных 
значений и от расчетных для тормозного момента (по математической модели). Вычислили процент ошибок, 
сравнив экспериментальный и математический результат. 

В таблице 6 приведены результаты подтверждения адекватности математической модели. 

















Таблица 6 
Результаты проверки адекватности модели 
Вязкость Плотность | Тормозной момент Тормозной момент Ошибки, 
№ Давление 

жидкости жидкости (расчет) (математическая модель) % 

1 1,85 1,15 38 490 495 1 
1,60 1,10 13 406 434 

3 1,45 1,12 26 423 447 5 


























Итак, результат подтверждает адекватность модели. Во-первых, уровень ее ошибочности оказался 
приемлемым. Во-вторых, выявленные ошибки имеют научное объяснение. Дело в том, что параметр плотности 
не зависит от вязкости, и это учтено в модели. Но эксперименты, проведенные ранее для подтверждения этой 
точки зрения, выявили между этими показателями корреляцию на уровне 6 % — и это примерно эквивалентно 
проценту ошибок, допущенных при использовании математической модели (в сравнении с расчетными 
данными). 

Для выражения совокупного влияния плотности и вязкости жидкости на работу гидравлического 
тормоза можно ввести в модель параметр и’р. Однако это необязательно, учитывая сравнительно низкий 
процент ошибок. В противном случае предложенная модель значительно усложнится. 

Имитация полевых условий. Предлагаемая математическая модель выражает физическое явление, 
выявленное и исследованное в лабораторных условиях при начальном тормозном моменте 355 Н\м. 

В полевых условиях балласты характеризуются начальным значением тормозного момента Мво. 
Предлагается ввести в математическую модель параметр, который не изменит ее форму, но, предположительно, 


отразит «полевой» балласт: 
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М, =58,5Р- 284, бр +1, 29 +10,62+ (Ми М»), (19) 
где Мв — тормозной момент «полевого» балласта после применения условий (давление, плотность и вязкость); 
Мв, — начальный тормозной момент для полевого балласта; Мз„ — начальный тормозной момент 


лабораторного балласта, равный 355 Н\м. 

Это основная гипотеза. Она обосновывается следующим образом. Размеры лабораторного ингибитора 
выбраны не случайно, а так, чтобы соответствовать реальным меньшим размерам тормозов, производимых 
компанией Рагтас [..[.С (модель 112—500)". 

Обсуждение и заключения 

1. Математическая модель создана с помощью метода обратного анализа, который связывает 
тормозной момент гидродинамического тормоза с рабочими параметрами (давление, вязкость и плотность) 
жидкости внутри балласта. 

2. Эксперименты при различных случайных значениях переменных доказали адекватность модели. 
Сравнивались значения тормозного момента, определенные экспериментально и с помощью созданной модели. 
Ошибочность не превысила 6 % (таблица 6). 

3. Предложена поправка к модели для симуляции работы гидродинамического тормоза в полевых 
условиях. 

4. Опираясь на полученные результаты, в будущих исследованиях целесообразно испытать созданную 
модель в полевых условиях с реальной полезной нагрузкой. 
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